Equilibrium magnetization at the boundary of a magnetoelectric antiferromagnet 
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Symmetry arguments are used to show that a boundary of a magnetoelectric antiferromagnet has 
an equilibrium magnetization. This magnetization is coupled to the bulk antiferromagnetic order 
parameter and can be switched along with it by a combination of E and B fields. As a result, the 
antiferromagnetic domain state of a magnetoelectric can be used as a non-volatile switchable state 
variable in nanoelectronic device applications. Mechanisms affecting the boundary magnetization 
and its temperature dependence are classified. The boundary magnetization can be especially large 
if the boundary breaks the equivalence of the antiferromagnetic sublattices. 



Magnetoelectric antiferromagnets (AFM) develop a 
magnetization M (or electric polarization P) in the bulk 
when an electric (or magnetic) field is applied [J-Q . This 
property is due to the presence of a magnetoelectric term 
in the free energy, Fme = ~OiikEiHi;, where is the 
magnetoelectric tensor; the latter is odd under time re- 
versal. An AFM is magnetoelectric if the presence of an 
invariant polar vector E can reduce its magnetic point 
group to a ferromagnetic one [1, 0] • 

Magnetoelectric and multiferroic materials can pro- 
vide the necessary response to allow electrical switch- 
ing of the magnetic state 0, [3-01 and potentially en- 
able fast, high-density, low-power, and non- volatile mem- 
ory devices (magnetoelectric memory) [8-ll|. To en- 
able easy readout of the magnetic state, the magneto- 
electric or multiferroic layer needs to be coupled to a 
proximate ferromagnetic layer. This coupling requires an 
exchange bias 12h14|| at the interface, which is the time- 
reversal-breaking shift the hysteresis loop of the ferro- 
magnet along the magnetic field axis. Much attention in 
this context has been focused on the room-temperature 
multiferroic BiFeOa, but the destabilizing effects of fer- 
roelastic strains and depolarizing fields need to be cir- 
cumvented for non- volatile operation [l5l |. Ferroelastic 
strain could be avoided by using a multiferroic material 
with linear coupling of P and M, but suitable materials 
for room-temperature operation are not available [3l ■ 

An alternative approach to electric magnetization con- 
trol uses the AFM order parameter of a magnetoelectric 
material as the switchable state variable. Magnetoelec- 



tric switching of Cr203 was shown 17[ to induce a re- 
versible switching of the exchange bias polarity in the 
proximate ferromagnetic Pd/Co multilayer on the macro- 
scopic scale. It was argued |17| that this effect is a man- 
ifestation of the equilibrium boundary magnetization of 
a magnetoelectric, which is required by symmetry and 
couples to the bulk AFM order parameter. Essentially, 
the boundary reduces the symmetry in a similar way to 
the electric field. Another manifestation of this effect is 
the spin polarization of the photoelectron current emit- 
ted from the free Cr203 (0001) surface fvi\ . 

Macroscopic signatures of boundary magnetization of 



Cr203 [17] show that the lack of macroscopic time- 
reversal symmetry in a magnetoelectric can translate into 
strong spin asymmetry at its boundary. However, the mi- 
croscopic mechanisms of this effect are not understood. 
In this Letter the salient features of boundary magnetism 
of magnetoelectrics are analyzed from the general point 
of view. A rigorous microscopic proof of the existence 
of equilibrium boundary magnetization is given, and its 
microscopic mechanisms are classified. In particular, it is 
shown that the effects can be very large if the boundary 
breaks the equivalence of the AFM sublattices. 

Consider a macroscopically flat boundary (surface or 
interface) of an AFM with an external normal n, which is 
allowed to have roughness and all possible terminations 
distributed with equilibrium Gibbs probabilities. The 
magnetic structure of the boundary is also assumed to 
be equilibrium, subject to the constraint that the bulk of 
the crystal is in the single AFM domain state [l8l |. 

We are generally interested in the response of the 
boundary on the macroscopic scale to an external probe 
which couples to the magnetic moment density m(r). 
This probe can represent spin-resolved photoemission, 
magneto-optic Kerr effect, exchange bias with a ferro- 
magnet, or magnetometry. For typical probes the mea- 
sured quantity is an odd functional G{mi{r)}, such that 
G{mi{r)} = G{mi{r + 1)} for any shift t (or at least for 
any translation vector of the bulk lattice treated as non- 
magnetic, such that t • n is not large compared to the 
equilibrium roughness) . Hereafter a probe is assumed to 
satisfy this condition [l9| . The component nii is selected 
by the polarization of the probe. 

The following arguments do not depend on the nature 
of the probe. For definiteness, let us select the boundary 
moment A4 as the probe, defined in a way that satisfies 
the above requirement of translational invariance. Specif- 
ically, if the magnetic unit cell is larger than the struc- 
tural unit cell, the magnetic moment A4 of the boundary 
region must be averaged over the different ways of sep- 
arating this region from the bulk, related to each other 
by purely structural translations (see Fig. 1). Surface- 
sensitive probes like exchange bias or spin-polarized spec- 
troscopies are free from this complication. 
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The macroscopic boundary magnetization is given by 
the equihbrium Gibbs average Mf, = (A^}/A, where A 
is the boundary area, and the thermodynamic hmit of 
large A is assumed. Its i-th component vanishes only if 
any boundary configuration (termination and magnetic 
structure) has an energetically degenerate one with a re- 
versed Mi] otherwise it is finite. Such degeneracy occurs 
only if the bulk magnetic space group of the crystal con- 
tains an operation under which the vector n is invariant, 
and Mi is odd. All degenerate boundary configurations 
can be generated by bulk space group operations; this is 
the boundary generalization of the Aizu procedure [20| . 
In particular, energetically degenerate atomic steps are 
automatically accounted for, as shown in Fig. 1 and 2b. 
However, since both n and Mi are invariant with respect 
to any translation, the latter can be disregarded, and we 
are led to consider only the magnetic point group. The 
presence of an invariant polar vector n selects the same 
subgroup of the bulk magnetic point group as a homoge- 
neous E field. It follows that the boundary acquires the 
same magnetization components as the bulk with an ap- 
plied E field in the direction of n. Therefore, the bound- 
ary develops finite equilibrium magnetization only if the 
bulk is magnetoelectric. This conclusion is equally valid 
for a metallic AFM whose magnetic point group would 
make an insulator magnetoelectric. 
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FIG. 1: Example of a bulk space group operation producing 
a degenerate configuration at the surface. An antitranslation 
TxR, where Tis a pure translation normal to the surface (red 
arrow) and R is time reversal, "cuts off" an atomic layer at 
the surface of a tetragonal lattice with layered AFM ordering. 
This bulk space group operation forbids both magnetoelec- 
tricity and equilibrium boundary magnetization. Black dot- 
ted line shows the boundary. Red and blue dotted lines show 
two inequivalent types of boundary/bulk partitioning, which 
are related by a non-magnetic translation T. The boundary 
moment averaged over these two partition types vanishes, be- 
cause the two types of atomic steps are degenerate. 



This conclusion reflects the fact that the free energy of 
the system with a boundary depends on the polar vector 
n as a macroscopic parameter. Just as in the bulk, the 
existence of the magnetization at the boundary can be 
deduced from the reduction of the bulk magnetic point 



group by the presence of the boundary, because transla- 
tions do not affect n or Mf,. From a different angle, a 
boundary can generate a magnetization only if its zero 
value is not protected by macroscopic time-reversal sym- 
metry in the bulk; this singles out the magnetoelectrics. 

Thus, equilibrium boundary magnetization of a mag- 
netoelectric is finite unless azk = for all k in the ref- 
erence frame where n is parallel to the z axis. If this 
magnetization is finite for the given n, it is also necessar- 
ily finite for any particular termination with this orienta- 
tion, because the magnetic symmetry group of the latter 
is a subgroup of the former. 

Probes that are not surface-sensitive, such as magne- 
tometry with G{mi{r)} = J mi{r)d^r, measure the sum 
of contributions of two film boundaries. The total mag- 
netization is non-zero unless there is a bulk symmetry 
operation interchanging the boundaries. It is always non- 
zero if these boundaries are with different materials. 

The exchange bias induced in a proximate ferromag- 
netic film by a magnetoelectric is fundamentally different 
from the conventional exchange bias, which is due to a 
small excess magnetic moment "frozen-in" in the AFM 
during field-cooling. In particular, this non-equilibrium 
character typically leads to an irreversible decline of the 
exchange bias as the magnetization of the ferromagnetic 
layer is repeatedly reversed — the so-called training ef- 



fect [I3|,ll|. By contrast, the switchable exchange bias 



observed in Ref. [17| is an equilibrium property and does 
not exhibit the training effect. 

Since the effect of the boundary is comparable to that 
of E of the crystal-field scale, even simple extrapola- 
tion suggests that the induced magnetic moments at the 
boundary can be a few orders of magnitude larger than 
those achievable due to the bulk magnetoelectric effect. 
In fact, some mechanisms do not contain any small pa- 
rameters and are capable of producing magnetizations of 
the order of a few Bohr magnetons per boundary site. I 
now classify these mechanisms. 

All mechanisms producing linear magnetoelectric re- 
sponse in the bulk 0] can generate boundary magneti- 
zation as well; these include changes in (A) the g-tensor 
(here we can also include hybridization-induced changes 
of the spin moment), (B) the single- ion anisotropy ten- 
sor, (C) the intrasublattice symmetric coupling (including 
Heisenberg exchange and dipolar interactions), and (D) 
the Dzyaloshinsky-Moriya [U exchange coupling in- 
duced by E or by the boundary. All of these except C 
involve relativistic terms in the Hamiltonian. 

Consider the usual case of a colhnear AFM. Mechanism 
C is active only if the perturbation breaks the equivalence 
of the AFM sublattices. This symmetry breaking can 
be identified by analyzing the so-called black-and-white 
(Heesch-Shubnikov) point group based on the decoration 
of the magnetic sites with Ising spin variables instead 
of axial vectors. If the perturbation (polar vector E or 
n) removes all symmetry operations mapping black and 
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white sites onto each other, the AFM sublattices become 
inequivalent; otherwise they do not. If the black-and- 
white point group allows Ising ferromagnetism, the true 
magnetic point group is also ferromagnetic, but the re- 
verse is not necessarily true. In the first case mechanism 
C is allowed, but in the second case the magnetoelec- 
tric response occurs only through spin canting due to a 
relativistic perturbation. Thus, certain components of 
the magnetoelectric tensor, and likewise the boundary 
magnetization for certain directions of n, may have no 
contribution from mechanism C. For example, consider 
Cr203 (magnetic point group 3m) with E or n oriented 
parallel to one of the three U2 axes or to one of the three 
ad planes bisecting them. The corresponding symmetry 
operation is not removed by E or n; since both U2 and 
ad interchange the AFM sublattices, the latter remain 
equivalent. However, the appearance of M parallel to E 
or Mb parallel to n through spin canting is allowed. 

If the equivalence of the AFM sublattices is broken by 
the boundary, the consequences are far more drastic than 
in the bulk mechanism C. For any particular boundary 
termination, even without reconstruction, the sites cor- 
responding to different AFM sublattices are structurally 
different. For example, all the sites closest to the bound- 
ary can have spins "down," while there is no equivalent 
termination with boundary spins "up." The situation 
is illustrated in Fig. 2b using Fe2Te06 (magnetic point 
group A/ mmm ) as an example. In this figure, termina- 
tions A (with boundary spins up) and B (with boundary 
spins down) are structurally distinct, and therefore they 
occur with different probabilities in equilibrium. Even 
if they did appear with equal weights, they are inequiv- 
alent magnetically, and don't generally add up to zero. 
Indeed, there are several mechanisms leading to the devi- 
ation of the average magnetic moments on the boundary 
sites from the bulk ones (see Fig. 2c-2e): (SI) Different lo- 
cal environment of the magnetic sites near the boundary 
leads to a different local magnetic moment, and perhaps 
even a different atomic multiplet. (S2) Since the trans- 
lational symmetry is broken by the boundary, any sym- 
metric exchange interaction (even purely intersublattice 
one) leads to different thermal averages at the boundary. 
(S3) The exchange coupling near the boundary can be so 
different from the bulk that the AFM ordering pattern 
can change to ferrimagnetic there. Mechanism SI can 
be viewed as a boundary analog of bulk mechanism A, 
and S2 is the boundary counterpart of mechanism C. SI 
and S2 are always present if the black-white symmetry 
is broken. Apart from these mechanisms affecting the 
magnitudes of the local moments, the coupling of these 
moments to the external probe can be different. For ex- 
ample, in the exchange bias setup, the sites closest to the 
boundary are expected to have the strongest exchange 
coupling to the proximate ferromagnet. 

The bulk linear magnetoelectric effect appears in the 
free energy expansion as a second-rank tensor aik- This 



is not the case for the equilibrium boundary magneti- 
zation. Since the free energy is a non-analytic function 
of n 2JI , the boundary magnetization, given by its field 
derivative, is also non-analytic. Thus, even if the AFM 
sublattices are equivalent for some symmetric directions 
of n, forbidding mechanisms S1-S3 for these orientations, 
these mechanisms can still generate large boundary mag- 
netization for less symmetric orientations. 

Boundary magnetization M{, vanishes at the bulk Neel 
temperature, but different mechanisms can be partially 
distinguished experimentally based on its temperature 
dependence. The thermal mechanism S2 can lead to a 
non-monotonic contribution with a maximum, similar to 
the bulk behavior of mechanism C (cf. azz in Cr203). 
Other mechanisms should lead to monotonically de- 
creasing with T. All these mechanisms do not contain 
any small parameters and can produce M;, of the or- 
der of a few Bohr magnetons per boundary site. The 
non-monotonic temperature dependence of the exchange 
bias field observed in the heterostructure of Ref. [l3l sug- 
gests that mechanism S2 plays an important role at the 
Cr2O3(0001)/Pd interface. 

Roughness-insensitive mechanisms based on 
Dzyaloshinskii-Moriya interaction at the compensated 
interface were proposed [2^ to explain the exchange 
bias induced by single-domain multiferroic BiFeOs in a 
proximate FM [lo| . They are, however, fundamentally 
different from all the mechanisms discussed here, because 
BiFeOa is weakly ferromagnetic in the bulk, and also 
its M and P are not linearly coupled [2^. Therefore, 
the boundary does not induce additional magnetization 
components that are forbidden in the bulk. 

Much attention is devoted to the search of a room- 
temperature multiferroic material with linear coupling 
between the electric polarization P and magnetization 
M, because its M could be switched along with P by elec- 
tric field only 0-0, S 13 ■ The paraelectric phase of such 
a material is a magnetoelectric AFM 3, 23] . The desir- 
able geometry involves voltage applied across the multi- 
ferroic film, with E normal to its surface Q . Equilibrium 
magnetization Mf, necessarily exists at the boundary of 
such a material. This Mf, has the same components as 
the bulk M coupled to P, but Mf, is coupled to the bulk 
AFM order parameter. Ferroelectric switching directly 
switches only the part of M linearly coupled to P, but 
not M;,. In addition, the states with parallel or antipar- 
allel Mf, and M are non-degenerate, meaning that one of 
them is metastable or even unstable. Thus, the presence 
of intrinsic boundary magnetization may hamper purely 
electric magnetization switching using a multiferroic with 
linear coupling of P and M. 

Equilibrium boundary magnetization of magnetoelec- 
tric AFMs has far-reaching consequences for the de- 
sign of magnetic nanostructures. First, very large ex- 
change bias fields should be achievable in magnetoelec- 
tric/ferromagnet bilayers, comparable to the estimates 
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FIG. 2: (a): Projected unit cell of the trirutile lattice assumed by magnetoelectric Fe2Te06 (schematic) pj]. Black circles: Te 
atoms; circles with arrows: Fe atoms with spin directions. O atoms, whose positions make the lattice non-symmorphic, are 
not shown, (b): Boundary termination cutting between the Fe layers (indicated as type A), including an atomic step. The 
same AFM sublattice appears at the surface (blue circles). The red arrow shows the bulk symmetry operation (4- fold screw 
rotation) mapping degenerate atomic steps onto each other. Termination type B (blue dotted line) puts another sublattice at 
the surface, but it is not related by symmetry and not degenerate with type A. (c)-(e): Mechanisms affecting the boundary 
magnetization. Termination type B is used as example; Te sites shown in gray may or may not be occupied, (c) Mechanism 
SI: Changed spin state (blue arrow), (d) Mechanism S2: Unequal thermal averages. Shades of blue and arrow size indicate 
the degree of thermal disorder, (e) Mechanism S3: Flipped spin direction (blue arrow). 



for a fully uncompensated AFM interface [12, 28|. Sec- 
ond, magnetoelectric AFMs are precisely the materials 
that can be switched between the time-reversed AFM do- 
main variants by a simultaneous application of E and B 
fields [l, , thereby switching the boundary magnetiza- 
tion and the exchange bias field . The B field may be 
permanent, while E may be provided by a voltage pulse 
across the magnetoelectric film. Since no depolarizing 
fields or elastic strains are involved, the AFM domains 
are stable, and this switching is fully non-volatile. Some 
device architectures based on the magnetoelectric active 
layer, where the operation is based on the linear bulk 
magnetoelectric response, were described by Binek and 
Doudin [§| . These architectures are greatly facilitated by 
the existence of a switchable equilibrium boundary mag- 
netization, which moreover allows the AFM domain state 
to be used as a switchable state variable. 

In summary, symmetry requires that magnetoelectric 
antiferromagnets possess a finite boundary magnetiza- 
tion in thermodynamic equilibrium. This magnetization 
is particularly large if the boundary breaks the equiv- 
alence of the antiferromagnetic sublattices; specific mi- 
croscopic mechanisms have been classified. This un- 
derstanding will hopefully stimulate further studies of 
boundary magnetization of magnetoelectrics and its ex- 
ploitation in nanoelectronic devices. 

This work was supported by NSF MRSEC, the Nano- 
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Corporation. 
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